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Abstract 

We have observed an environmentally induced quantum dynamical phase transition in the dynamics of a two spin 
experimental swapping gate [J. Chem Phys. 124, 194507 (2006)]. There, the exchange of the coupled states ||, J.) 
and ||, T) gives an oscillation with a Rabi frecuency b/h (the spins coupling). The interaction, H/tse with a spin-bath 
degrades the oscillation with a characteristic decoherence time. We showed that the swapping regime is restricted 
only to 6rsB > h. However, beyond a critical interaction with the environment the swapping freezes and the system 
enters to a Quantum Zeno dynamical phase where relaxation decreases as coupling with the environment increases. 
Here, we solve the quantum dynamics of a two spin system coupled to a spin-bath within a Liouville-von Neumann 
quantum master equation and we compare the results with our previous work within the Keldysh formalism. Then, 
we extend the model to a three interacting spin system where only one is coupled to the environment. Beyond a 
critical interaction the two spins not coupled to the environment oscillate with the bare Rabi frequency and relax 
more slowly. This effect is more pronounced when the anisotropy of the system-environment interaction goes from 
a purely XY to an Ising interaction form. 
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The control of quantum dynamics started to receive 
much attention because of its relevance to several appli- 
cations ranging from quantum information processing 
[1] to nanotechnology [2]. Since the environment usu- 
ally acts degrading the quantum dynamics of a system 
[3], its decoupling becomes a major challenge. Many 
techniques [4,2,5] are developed to avoid this loss of in- 
formation that is characterized by a decoherence rate 
l/r^. Spin systems are ideal candidates to test the pro- 
cedures for QIP. Recently, we observed an environmen- 
tally induced quantum dynamical phase transition in 
the dynamics of a two spin experimental swapping gate 
[6]. There, the exchange of the coupled states |t,4) 
and ||, T) gives an oscillation with a Rabi frequency 
b/h (the spins coupling). We showed that the swapping 
regime is restricted only to brsE > h, where H/tse is 
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the system-environment (SE) interaction, and essen- 
tially within this regime 1/t0 ~ 1/tse- However, be- 
yond a critical interaction with the environment, the 
swapping freezes and the system enters to a Quantum 
Zeno dynamical phase where the relaxation rate de- 
creases as the coupling with the environment increases 
l/ 1 "^ ~ & 2 T"sE- Here, we will show how this criticality is 
useful to "isolate" a spin-pair. 

Firstly, we solve the quantum evolution of a two spin 
system coupled to a spin-bath within the Liouville-von 
Neumann quantum master equation [7,8] , and compare 
the result with our previous one [6] using the Keldysh 
formalism [9]. Then, we extend the model to a three 
interacting spin system where only one is coupled to 
the environment. Beyond a critical interaction, the two 
spins not directly coupled to the environment oscillate 
with their bare frequency and relax more slowly. In a 
two spin system there is always a critical point that de- 
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pends on the anisotropy relation of the SE interaction 
quantified as the ratio between Ising and XY terms. 
However, in the three spin system, the decoherence rate 
has a smooth cross-over from proportional to the SE 
interaction to inversely proportional to it. This cross- 
over approaches a critical transition as the anisotropy 
of the SE interaction goes from purely XY to Ising. 

Experiments in Ref. [6] focus on two interacting spins 
1/2 coupled to a spin-bath as modeled by the Hamil- 
tonian H = Hs + Hse + He, with 

Hs = ho;* {II + IS) + b (7+/ 2 - + /r/ 2 +) /2, (1) 
Use = aI%F z + (3 (l+F~ + I^F+) /2, (2) 

He = !^>^^2<Ji 

+ E 2<i <^ [2/f /; - 1 (/+/- + it it)} , (3) 

where Hs is the system Hamiltonian of the two cou- 
pled spins, He is the spin-bath Hamiltonian with a 
truncated dipolar interaction and Hse is the SE in- 
teraction with F u = 53 j b 2 il", u — x,y,z and F = 
(F x ± iF v ) . Hse is an Ising interaction if /3/a = and 
a XY, isotropic (Heisenberg) or the truncated dipolar 
interaction if a//3 = 0, 1, —2 respectively. We use the 
model proposed by Miiller et al. [10] to calculate the 
spin dynamics in the experimental system of Ref. [6] . 
The model assumes that only one spin interact with 
the spin-bath which is described in a phenomenological 
way. In a quantum mechanical relaxation theory the 
terms F w are bath operators while in the semi-classical 
theory [7,8] F u (t) represent classical stochastic forces. 
As the experimental conditions justify a high tempera- 
ture approximation, the semiclassical theory coincides 
with a quantum treatment. By tracing on the bath 
variables the random SE interaction Hamiltonian is 

Hse (t) = aF z (t) J| + (3 [F~ (t) 1+ + F + (t) / 2 "] . 

(4) 

The time average of these random processes satisfy 
F u (t) = 0, where their correlation functions are 
g (u,v) ( T ) = pu ^ pv , (t + T y Following the usual 

treatment to second order approximation, the dynam- 
ics of the reduced density operator is given by [7,8] 
d 
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[Hs,<T{t)]-T{<r(t)-<T }, (5) 



where the relaxation superoperator T contains the SE 
interaction. It accounts for the dissipative interactions 
between the reduced spin system and the spin-bath and 
imposes the relaxation of the density operator towards 
its equilibrium value er . We assume that the correla- 
tion times of the fluctuations are extremely short com- 
pared with all the relevant transition rates between 
eigenstates of the Hamiltonian. In this extreme narrow- 
ing regime or fast fluctuation approximation we obtain 

r W = l*Etu.«J {u ' v) (o) [is, [i v 2 ,*]] , (6) 



where J (u ' v) (w) = f™ x drg^ v) (r) exp {-icur} is the 
spectral density and = (aS u ,z + (3 (S u ,x + $u,y)) 
(a5 V;Z + (3 (5 V;X + 5 Vz y)). The spatial directions are sta- 
tistically independent, i.e. g( u ' v ^ (r) = if u 7^ v. No- 
tice that the axial symmetry of Hs around the z axis 
leads to the impossibility to evaluate separately J x and 
J y , where J u = \ j {u ' u) (0). Thus, they will appear 
only as the averaged value J xy = {J x + J y ) /2. The 
superoperator now can be written as 

r{a} = r zz [I2AI2M] 

+ r XY ([If,[/f, a]] + [/!,[/!, a]]), (7) 

where Fzz = ct 2 J z and Txy = ft 2 Jxy Note that Tzz 
and Txy contain the different sources of anisotropy. 
The usual approximation considers J x — J y — J z 
(identical correlations in all the spatial directions) 
and a = (3 = 1 (isotropic interaction Hamiltonian) 
[10]. A better approximation considers a dipolar in- 
teraction Hamiltonian, i.e. a = —2/3 = 2 [11,6]. We 
consider the experimental initial local polarization on 
site 2, er(0) = [1 + /fefiwoj-TI] /Tr {1} and the spin- 
bath polarized, where /3b = 1/ (ksT). As the final 
state reaches the temperature of the spin-bath, ao = 
[1 + (3-g.huJoi (If + Jf )] /Tr {1}. Here, ao commutes 
with Hs, not containing coherences with AM > 1. 
Following the standard formalism [7,8], we write the 

superoperator T using the basis of eigenstates of the 
system Hamiltonian (1). After neglecting the rapidly 
oscillating non-secular terms with respect to the Hamil- 
tonian, i.e.: Tzz,rxY <§C b, we solve Eq. (5) and we 
calculate the magnetization of the spin 1 obtaining an 
extension of the result of Ref. [10] . Our essential contri- 
bution is that we specifically account for the anisotropy 
arising from the nature of SE interaction reflecting it 
in Tzz = a 2 Jz and r XY = f3 2 J xy , 



M n (t) = Tr{I z l a{t)} = M 



i-i.- 



-icosHt) e -( 2 " 2 ^+« 2 ^)'/ 2 



, (8) 



where loo = b/h and Mo = [3&TUi>oi/&. 

If we relax the condition F Z z, Txy <S b, i.e., we do 
not neglect the non secular terms for the superoperator 

T, the dynamics still occurs in the Liouville space of 
the populations and ZQT. We obtain exactly the same 
solution derived in [6] within the Keldysh formalism: 

Mq (t) 



Mo 



(l - aoe- Rat - 01 cos [(w + LR 2 ) t + (j>o] e _fllt ) 

(9) 



where the real functions to, Ro, Ri and R2 as well as ao, 
ai and <f>o depend exclusively on b, 1/tse = T Z z + T X y 
and pxy = Txy/ (rzz + Txy). The complete analyti- 
cal expression is given in Ref. [6] . There, we showed that 
if the SE interaction is anisotropic (Fxy =fc Tzz), the 
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functional dependence of to on tse and b yields a criti- 
cal value for their product, brss/h = k PXY , where the 
dynamical regime changes. We called this phenomenon 
a quantum dynamical phase transition [6,9] ensured 
by the condition L0R2 = 0. This is evidenced by the 
functional change (non-analyticity) of the dependence 
of the observables (e.g. the swapping frequency to) on 
the control parameter brss/h. This non-analyticity is 
enabled by taking the thermodynamic limit of an infi- 
nite number of spins. [12]. One identifies two paramet- 
ric regimes: 1- The swapping phase, which is a form of 
sub-damped dynamics, when &TSE/fc PXY > h (R2 = 
in Eq. (9)). 2- A Zeno phase, with an over-damped dy- 
namics for frrsE/fcpxY < ^ arising on the strong cou- 
pling with the environment (zero frequency, i.e. lu = 0, 
in Eq. (9)). In the neighborhood of the critical point 
the swapping frequency takes the form: 

_ f s pxy \J ( b / h ) 2 - ^pxyAse ^SE/fcp XY > h 

\ brsE/fcpxY ^ ^ 

(10) 

The parameters a PXY an d k PXY only depend on pxy 
which is determined by the anisotropy of the interac- 
tion Hamiltonian. 

Three spin system. The system Hamiltonian is 

Hs = Huj l (/ z + /i z + /|) 

+ b (/+/r + j - /+) /2 + & (/+/ 2 - + ir/ 2 + ) /2, (11) 

and the environment and SE Hamiltonian remains as 
before. Also, the environment is coupled to only one 
spin of the system. We solve Eq. (5) as before, with- 
out neglecting non-secular terms of the relaxation su- 
peroperator. Considering the initial condition a (0) = 
(1 + f3hcooilo) /Tr{l} and a polarized spin-bath, the 
magnetization at site is 

Mi S (t) = Tr{J ( ^a (t)} = M [l - a e~ Rot - aie~ Rv 

+a 2 sin (oj 2 t + fa) eT R ' 2t + 03 sin (w 3 t + fa) e~ fl ' 3t ] . 

(12) 

The coefficients Oj , Ri, u)i and fa are real and they are 
functions of 6, 1/tse andpxY- IfpxY / the final state 
has all the spins polarized because a net transfer of 
magnetization from the spin-bath is possible. However, 
for an Ising SE interaction, pxy = 0, we obtain that 
Ro — and 1 — ao = 1/3 (the asymptotic polarization) 
because the final state is the quasi-equilibrium of the 
3-spin system [11]. In Fig. 1 we show the frequencies u>2 
and and the different relaxation rates as a function 
of (6rsE/ft) _1 when the SE interaction is Ising (jjxy = 
0) . Two changes, resembling the critical behavior of two 
spin systems are observed. The same phenomenon oc- 
curs in Figure 2 a) where the coefficients o» are shown. 
The polarization evolution of an isolated 3-spin sys- 
tem is Miz (t) = Ma [3 + 4 cos + cos (wgi)] where 



uj? 2 = (\/2/4) b/h and = (y/2/2) b/h are the natural 
frequencies. When (bTss/h)' 1 << 1, we observe that 
lo 2 — » W3 — > W3, ai — > 1/3 — 3/8, 02 — » 1/2 and 
03 — > 1/8 as expected for an isolated 3-spin dynamics. 
The dependence of ^3 as a function of (fersE/ft) 1 is 
similar to that of the swapping frequency in Ref. [6]. 
However, instead of becoming zero when the SE inter- 
action increases, it suddenly stabilizes at uio = b/h, 
the bare 2-spin Rabi frequency. At the same point u)2, 
R2 and R3 also have a sudden change. While R 2 and 
R3 initially grow as (fcrsE/fi) -1 , there R2 increases the 
growing speed while R3 begins to decay as in the Zeno 
phase of Ref. [6]. Moreover, the coefficients 02 and 03 
switch between them, 02 suddenly goes down and ft3 
goes up. These coefficients are the weights of the differ- 
ent frequency contributions in the time evolution. The 
changes on the decoherence rates and on the weight co- 
efficients of the different oscillatory terms beyond the 
critical interaction (region) lead the system to oscillate 
with the bare Rabi frequency of two spins decoupled 
from the environment. If we keep increasing the con- 
trol parameter (brsE/h)~ l , this effect is enhanced by 
the next transition. After the second transition, R\ be- 
gins to decrease as in the Zeno phase behavior of Ref. 
[6]. As the term of Eq. 12 that relaxes with Ri leads 
the system to the 3-spin quasi-equilibrium, when Ri 
goes down, this final state is approached much slower. 
The effect is to avoid the interaction between the 2- 
spins not coupled to the environment and I2 ■ After the 
second transition, the coefficient a\ goes abruptly to 
zero leading to a more pronounced "isolation" of the 2- 
spins. Thus, two dynamical regimes are observed: One 
characterized by the 3-spin dynamics for (brsE/h) _1 < 
1 and a second one, for (hrse/h)" 1 > 1, which has a 
2-spin behavior. 

Fig. 2 b) and c) show the temporal evolution of the 
magnetization of Eq. (12) on the 3-spins and the 2- 
spin regimes respectively. While in Fig. 2 b) the two 
frequency contributions are evident, in Fig. 2 c) only 
the bare Rabi frequency is manifested. In each graph 
we show two curves with different SE interactions. In 
Fig. 2 b), we show that increasing the SE interaction 
the decoherence rates increase. However, in the 2-spin 
regime (Fig. 2 c)) when the SE interaction is increased, 
the decoherence rate decreases leading to a better "iso- 
lation" . It is important to take into account that while 
the relaxation rates go to zero smoothly the swapping 
frequency acquire the bare value near the critical point. 
Another fact to remark is that this effect is more pro- 
nounced when the anisotropy of the SE interaction is 
close to a pure Ising SE interaction while an increase 
in the XY nature leads to a further smoothing of the 
transition. The reason is that, when pxY =fc 0, there is 
a net transfer of magnetization to the system which is 
redistributed between the three spins, this redistribu- 
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Fig. 1. (color online) a) Frequencies involved in the time evo- 
lution of the polarization in the 3-spin system as a function of 
(brsE/K) ~ 1 . Dashed lines represent the isolated system. Dot 
line corresponds to two spins decoupled from the environment, 
b) Different relaxation rates of the polarization. 




Fig. 2. (color online) a) Coefficients (weights) of the different 
terms of Eq. (12). At the critical region there is a switch 
between the 2-spin and the 3-spin regime, b) and c) Temporal 
evolutions of the polarization in the 2-spin and 3-spin regimes 
respectively for different tse- In b) b/h — 2tt X 1 kHz and 
tse — 1.43 ms for the thick line and tse — 10 ms for the thin 
line. In c) b/h — 2ir X 1 kHz and tse — 0.1ms for the thick 
line and tse — 0.01 ms for the thin line. 



tion begins to be slower at the second transition when 
i?3 goes down. In contrast, for a purely Ising interac- 
tion, there is no net polarization transfer and a purely 
decoherent process at site 3 freezes its dynamics but its 
fast energy fluctuations prevent the interaction with 
the other spins. 

In summary, we found an analytical expression for 
the two spin dynamics plus a spin-bath of the experi- 
mental swapping gate [6] , showing that standard den- 
sity matrix formalism leads to a quantum dynamical 
phase transition as does the Keldysh formalism [6,9]. 
Here, we extended the model to a 3-spin system and 
showed that beyond a critical region the two spins be- 
come almost decoupled from the environment oscillat- 
ing with the bare Rabi frequency and relaxing more 
slowly. While in the two spin swapping gate the dy- 
namical transition is critical, when we extend the sys- 
tem to 3-spin the criticality is smoothed out. How- 
ever, enough abruptness remains to give the possibility 
to use it to "isolate" a two spin system with a finite 
system-environment interaction. 
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